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Abstract. When combining the observations of spatial 
abundance and two-point correlation function of QSOs, 
we can effectively set constraints on models of cosmic 
structure formation. Both the abundance of gravitation- 
ally confined halos and their two-point correlation func- 
tions can be calculated in the conventional Press-Schechtcr 
formalism. We apply this method to examine the proper- 
ties of possible host halos of QSOs in three popular mod- 
els: the standard cold dark matter (SCDM) model, the 
low density flat cold dark matter (LCDM) model and the 
cold-plus-hot dark matter (CHDM) model. The LCDM 
and CHDM models are normalized to the COBE-DMR 
observations, and the SCDM is normalized to cr = 0.58. 
We find that the SCDM and LCDM models can pass the 
abundance-plus-correlation test for QSOs. However, the 
CHDM are difhcult to produce host halos to fit with the 
number density of high rcdshift QSOs and their clustering 
on large scales (10 h^^ Mpc) simultaneously. We stud- 
ied various mechanisms, originated both gravitationally 
and non-gravitationally, which may lead to a biasing of 
the halo clustering. We conclude that these effects are too 
weak in order to release the trouble of the CHDM models. 

Key words: quasar: general - cosmology: large-scale 
structure of universe - dark matter 



1. Introduction 

Recently, we have developed a method to test popu- 
lar dark-matter models of structure formation using the 
number density and the two-point correlation function of 
high redshift objects (Bi & Fang 1996). We found an ap- 
proximate expression of two-point correlation function of 
mass and collapsed halos in the Press-Schechter formalism 
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(Press & Schechter 1974, hereafter PS). In this approxima- 
tion, the nonlinear gravitational interaction was treated as 
the sum of various individual spherical top-hat clustering. 
These top-hat spheres consist of both collapsed PS halos 
and uncoUapsed regions. Moreover, the bias that massive 
PS halos have stronger correlation than the background 
mass can naturally be introduced by considering that no 
collapsed halo of mass M exists in initial regions (or top- 
hat spheres) of mass less than M. 

This method was applied to CIV absorption systems in 
QSO spectra. Because CIV systems should be hosted by 
collapsed halos, one can obtain an lower limit to the spa- 
tial number density of these host halos from the observed 
CIV number density. This requires that the hosts of CIV 
should consist of halos with mass as low as Mth- On the 
other hand, the two-point correlation of halos is stronger 
when the mass of the halos is larger. The observed two- 
point correlation functions set a lower limit, Mco, to the 
mass of host halos. Obviously, a reasonable model has to 
give Mth > Mco- The standard cold dark matter (SCDM) 
model normalized to tr = 0.58 and the low density flat 
cold dark matter (LCDM) model can pass this test. How- 
ever, for the cold-plus-hot dark matter (CHDM) models 
with parameters = 0-7 and flf, = 0.3, or fie = 0.8 
and flh — 0.2, the two-point correlation functions of ha- 
los with mass Mth are too small to explain the observed 
correlation functions. In order to have enough number of 
collapsed halos to host CIV systems, Mth should not be 
larger than 10^^ Mq. But the observed two-point corre- 
lation function on the scales of Av ~ 300 — 1,000 km/s 
indicates that Mco should not be less than 10"'^^ Mq. 

In this paper, we study the same problem, but us- 
ing QSOs as the discriminator. The topic of finding con- 
straints on dark matter models from QSOs is not new. As 
early as 1980s, the highest redshift of QSOs was used to 
rule out the model of hot dark matter (HDM), because the 
HDM predicted that collapsed halos cannot form at red- 
shift larger than 3, which is much less than the redshifts of 
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many existing QSOs. Consequently, the formation of high 
redshift QSOs cannot be explained in the HDM model 
(e.g. Efstathiou & Rees 1988). This redshift test is passed 
for current dark matter models. For instance, the tilted 
CDM, scaled CDM and CHDM models all are able to 
produce with the observed abundance of QSOs (Nusser 
& Silk 1993). However, the QSO abundance test alone is 
not free from the uncertainties of discriminating among 
the models like SCDM, LCDM and CHDM. While some 
works claimed that the CHDM model is consistent with 
the QSO abundance, other works found oppositely (Ma & 
Bertschinger 1994). We will show that adding a cluster- 
ing test to the abundance fitting will reduce the number 
of free parameters in the discrimination, and gives much 
better results. 

1.1. Abundance of QSOs 

The spatial number density of collapsed halos can be cal- 
culated from the standard PS theory. We define 6{x, z) to 
be the 3-D density fluctuation field of dark matter extrap- 
olated to redshift z assuming linear evolution. A density 
field (5/{(x, z), representing the smoothed fluctuation on 
scale R, can be derived from i5(x, z) by 

5fl(x, 2) = ^ / ^(x'. zW{R; x' - x)dx', (1) 

where the function VF(i?;xi — x) is the top-hat window 
for the comoving volume Vr = AttR^ /2>. The total mass 
within Vr on average is M = VrPo, where po is the mean 
density at the present if the scaling factor of the universe 
is set to be unity at 2; = 0. 

For Gaussian perturbations, the mass fluctuation 
within a top-hat window of radius R is described by vari- 
ance (7^(i?, z), which is determined by the initial density 
spectrum P{k) and normalization factor cts, i.e. cr(i? = 8 
h~^Mpc). In the case of f7 = 1 Einstcin-de Sitter universe, 
the linear evolution of the variance is cr^(i?, z) oc {l-\-z)~'^. 
Thus, the fraction of the total mass po^x having fluctua- 
tions larger than a given 5c in an arbitrary spatial domain 
Ax is 

Fr= r -^=^ expf ^3_\d5R. (2) 

If we take Sc{z) to be the critical overdensity for the 
collapse of the spherical mass M at z, Fr- po Ax should be 
identified as the sum of masses of all collapsed halos, each 
of which is massive greater than M. For the Einstein-de 
Sitter universe, we have Sc{z) = (1 -1-^)1.69. For the flat A 
universe, function Sdz) is calculated by Bi & Fang (1996). 
The differential —-^{FRpQAx)dM gives the total mass of 
collapsed halos in the mass range M to M + dM. Hence, 
if nc{M, z)dM is deflned as the spatial number density of 
halos between M and M + dM at z, we have 

-^iFRpoAyi)dM = nc{M, z)dM ■ MAx, (3) 



where we use the subscript c in ric to emphasize that it is 

for collapsed halos. Considering the cloud- in-cloud prob- 
lem, the above defined Uc should be multiplied by a factor 
of 2. The normalization Jg°° n^M, z)MdM = po can then 
be fulfilled. Therefore, the spatial number density of halos 
with mass between M and M + dM at redshift z is given 

by 

ne(M, z)dM = -^^edc ( -J^i^) dM, (4) 
' ^ MdM \V2a{R,z)J ^ ' 

where erfc(x) is the complementary error function. 

We win study models of SCDM, LCDM and CHDM, 
for which the density parameter f^o, the cosmological con- 
stant Ao, and erg) sje listed in Table 1 (the Hubble constant 
is taken to be Hq = 100 h km s~^ Mpc~^.) It is worth to 
point out that the as of the CHDM and LCDM models are 
compatible with the COBE-DMR observation, but that of 
the SCDM model is not. The value as = 0.58 gives a good 
fit for the SCDM model to all data except for the COBE 
result. We assumed that the initial primordial power spec- 
trum is the Harrison- Zel'dovich type. The linear transfer 
functions of the SCDM and the LCDM models are taken 
from Bardeen et al. (1986) and that of the CHDM model 
from Klypin et al. (1995). 

The cumulative number density, Nc{M), of all halos 
with mass greater than M should be 

/•CXD 

Nc{M,z)= / nc{M,,z)dM^. (5) 

JM 

The abimdance of halos calculated from Eq. (5) has been 
verified by a number of N-body simulations (e.g. Lacey & 
Cole 1994). 

The redshift evolution of the number density of col- 
lapsed halos, Nc{M,z), are shown in Figs, la, lb and 
Ic for the SCDM, LCDM and CHDM models, respec- 
tively. The eight curves in Figs, la and lb correspond 
to M = loi3-5+"0.25 ^^^Yi n = 0, 1...7 from top 

to bottom. In Fig. Ic, it is M = loii-5+"0.25 j^^^ .^j^j^ 
n = 0,1...7. 

Pairs and multiples of QSOs are not common to see. 
Therefore, it is reasonable to assume that each collapsed 
halo hosts only one, or at most a few QSOs. The number 
density of QSOs with magnitude Mb < —26 is plotted 
as crosses in Figure 1. The data points are taken from Pei 
(1995) whose result is based on the observations of Hewett 
et al. (1993) and Schmidt et al. (1992). In Pel's paper, n{z) 
is measured in the Einstein-de-Sitter cosmological model 
{Hq = 1, Ao = and h = 0.5). When comparing the ob- 
servations to the LCDM model, we have made corrections 
of the cosmological effect on n{z) due to the non-zero Aq. 

Fig. 1 shows that for all the models, the lower mass lim- 
its are mainly determined by the QSO abundance at z = 4. 
The possible hosts of QSOs in the SCDM and LCDM mod- 
els can be provided by halos with masses above lO^^ ''^, 
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Table 1. Parameters of the cosmological models 





SCDM 


LCDM 


CHDM 


h 


0.5 


0.75 


0.5 


no 


1 


0.3 


1 (fi^ = 0.3) 


Ao 





0.7 





(T8 


0.58 


1 


0.55 



and 10"'^^ Mq, respectively. In the CHDM model, the num- 
ber densities of collapsed halos at these mass range are 
much fewer. Therefore, for the CHDM, we have to include 
lower mass halos to host QSOs. The mass of possible QSO 
halos will be as low as about M — lO^^'^'"' Mq, which is 
more than two orders of magnitude below those in the 
models of SCDM and LCDM. 



a.SCDM 



-6 - 



-10 




2. Two-Point Correlation Functions of QSOs 

The spatial number density of uncollapsed spherical re- 
gions can also be calculated in PS formalism. We define 
m{Mo,r, z)dModr to be the number density of regions 
with mass Mq — > Mq + dMo and radius r ^ r + dr a,t 
z. Similar to Eq.(4), we have 



m{Mo,r,z) 



1 92 



Po 



Mo 2 dModr 



-erfc 



S{r,z) 



V2a{ro,z) 



(6) 



where ro is defined by Mq = Vopo = 47rro/>o/3. S{r,z) is 
the overdensity, for which a region Mq will evolve into a 
sphere with radius r at z. Function 6(r, z) is also calculated 
in Bi & Fang (1996). The total number of such spherical 
regions in an arbitrary volume Ax should on average be 
given by m(Mo, r, z)dModrAx. 

The mass correlation function ^(r) is determined by 

the relative enhancement of mass density in the spherical 
shell r ^ r + dr around area Ax. Only the spheres with 
radius r ^ r+dr can contribute to this enhancement. The 
mean enhancement of each Mq sphere is approximately 
described by its mass variance MQa'^{ro, z). Therefore, the 
mass correlation function can be estimated by 



/o m(Mo, r, z) dModr A:icM^a^{ro, z) 



f 

Jo 



poAnr'^dr ■ poAx 
dMom{Mo, r, z)M^a^{ro, ^)/47rrVo- 



(7) 



The factors poAnr^dr and poAx. are the mean mass in 
the spherical shell r — > r + dr and the volume Ax, re- 
spectively. In deriving Eq. (7) , we implicitly assumed that 
there is no correlation among the initial spheres, so the 
mass fluctuation in the shell is simply given by the sum 
of the individual components. 



b.LCDM 




Fig. 1. Evolutions of the comoving number density 
of halos with mass larger than a given number M. a.) 
for SCDM model, the eight curves correspond to M = 
2oi3.5+n0.25 ^^^^ ^-^^^i n = 0,1. ..7 from top to bottom, 
b.) The same as a.) but for LCDM model, c.) for CHDM 
model, here M = ioii-5+"0.25 ^j^j^ n = 0, 1...7. 

Because it is impossible that a collapsed halo with 
mass M formed in a uncollapsed region with Mq < M, 
the number density of M collapsed halos in an uncollapsed 
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sphere of mass Mo should be zero if M is greater than Mq. 

Thus, it is reasonable to assume that the number density 
of collapsed halos in each PS sphere r at z on average is 

M.(M,Mo)-|g forM>Mo, 

where the constant >1 > 1 is introduced to maintain the 
normalization condition 

rKVon{Mo, z)dMo = N,{M, z). (9) 
Jo 

We have then A = 2/eiic[S{R, z)/a{R, z)]. 

Like calculating Eq.(7), let's consider a typical spher- 
ical shell r ^ r + dr. Only the spheres with radius 
r ^ r + dr can contribute to the mass enhancement in 
this shell. The total number of micollapsed regions with 
radius r ^ r + dr and masses Mq Mq + dMg in a 
volume Ax is m(Afo, r, z)dModrAx. The total number of 
collapsed halos M in each Mq sphere is Afc{M)Vo, and the 
variance of the number is (A/'cVo)^(T^(ro, z). Therefore, the 
correlation function of collapsed halos with mass larger 
than M can be approximated as (Bi & Fang 1996) 

^{r;>M,z) = 

dMomjMQ, r, z)drAx ■ {J^^Vofa^ro, z) 
NAirr'^dr ■ NcAx 

/>00 

/ dMQm{Mo,r,z)V^a^{rQ,z)/A-Kr^ (10) 

JM 

Eq.(lO) implicitly assumed that in uncollapsed regions, 

the collapsed halos have the same linear variance as the 
mass. The approximation given by Eq.(lO) is found to be 
in good agreement with the linear approximation on scales 
larger than R, but higher than the empirical formalism of 
Hamilton et al. (1991) on scales less than R, here R is the 
scale within which the non- linear effects are significant. 
Therefore, Eq.(lO) can be used at least as a upper limit 
to the correlation function. 

Using Eq.(lO), we have calculated the correlation func- 
tion ^(r, > M) at redshift z = 2.5. The results are shown 
in Fig. 2, in which 2a and 2b are for the models of SCDM 
and LCDM, and 2c for CHDM. The eight curves in Figs. 
2a and 2b are corresponding to masses 10^'' '''+""'^'''AfQ, 
n = 0, ...7 from left to right, respectively. In Fig. 2c, it is 
10"-5+"°-25Mq and n = 0, ...7. 

The two-point correlation function of QSOs is found 
to be obey the same power law as galaxies ^gg(r) = 
AqqT^^'^ = {r/ro)~^'^, and the amplitude Agq 25 
(or the correlation length tq ^ 6 h~^Mpc) at z = 1.5 
when qo is taken to be 0.5 (Mo & Fang 1993). This gives 
^gg(r) > 0.1 for r = 5 - 10 h-^ Mpc. The clustering of 
QSOs on scales of r ~ 30 h~^Mpc is also foimd to be sig- 
nificant (Deng et al. 1994). Some observations and statis- 
tics have even indicated the possible existence of groups of 
QSOs with comoving sizes as large as about 100 h~^Mpc 



(e.g. Clowes & Camppusano 1991, Romberg et al. 1996.) 

Despite Aqq is decreasing for z > 1.5, the amplitude 
is still larger than 0.1 (Mo & Fang 1993; Romberg et al. 
1994.) 




5 10 15 

Comoving Distance r (ii"' Mpc) 




5 10 15 

Co,no>-,„s Dislaiice r (!,-' Mpc) 




5 10 15 

Comoving Distance r (1)-' Mpc) 



Fig. 2. Two point correlation functions of M halos at 
z = 2.5. a.) For SCDM model, the eight curves from left 
to right are corresponding to masses M = 1O^'^ ''^+"°'^^M0 
and n = 0, ...7, respectively, b.) The same as a.) but 
for LCDM model, c.) For CHDM model, here M = 
1O"-5+»o-25M0 and n = 0, ...7. 
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Fig. 2a and 2b show that the correlation function of 

M > lO^^'M© hales in the SCDM and LCDM models 
are consistent with the observational data. These halos 
have also proper spatial number density as the observed 
QSOs. The difference between the two models are very 
small. In fact, the top-hat evolutions of a spherical mass 
in SCDM and LCDM are indistinguishable, because it has 
almost the same dynamic trajectory in all flat universes 
with Ao < 0.8. 

On the other hand. Fig. 2c shows that the two-point 
correlation functions of the CHDM model arc very small 
on large scales. For halos with mass lO^^'^^M©, the am- 
plitude of the correlation function is well less than 0.1 
on scales of r =5 -10 h"^ Mpc. Therefore, the CHDM 
model seems to be unfavored under the abundance-plus- 
correlation test. According to the abundance, the mass 
of QSO hosts in the CHDM model should be as small 
as 1O^^M0. On the other hand, the two-point correlation 
function of such halos is much less than what is observed. 

3. Biasing problem 

A possible way to save the CHDM model is to assume that 
there is a very high biasing factor to raise the correlation 
function of 10^^-^Mq halos in the CHDM model. From 
Fig. 2c, if this factor is as large as 5, the power of two-point 
correlation functions mass would be able to fit 

with QSO's clustering. However, biasing is actually not a 
free parameter which can be arbitrarily chosen (Einasto 
et al. 1994). 

Generally speaking, the correlation functions of objects 
identified from a mass field will be different from the mass 
itself. A different method in the identification procedure 
can lead to a different biasing. The question now is: can 
we choose other QSO identifications to produce larger bi- 
asing, and then to give higher amplitude of the two-point 
correlation functions of QSOs in the model? 

As an example, let us first consider the QSO identifica- 
tion by velocity dispersion. It has been known for a decade 
that QSOs with low redshift are preferentially located in 
small groups of galaxies. This is in evidence from QSO- 
galaxy correlation function (Yee & Green 1987), CIV- 
associated absorption in high redshift radio-loud QSOs 
(Flotz et al. 1988), clustering analyses of QSO distribu- 
tion (Bahcall & Chokshi 1991) and the galaxy environ- 
ments around QSOs (Ellingson et al. 1991a). It has been 
also shown that the velocity dispersion of galaxies around 
QSOs is ~ 400 km s^^ (Elhngson et al. 1991b). There- 
fore the environment suitable of QSOs formation seems 
to be small groups of velocity dispersion ^ 400 km s"-'^. 
The strength of the QSO correlation function, is interme- 
diate between galaxies and rich clusters, so once again it 
is similar to galaxy groups. 

If we assume that high redshift QSOs formed in the 
same environment as that in the low redshift, the QSOs 
should be identified as collapsed halos with 3-dimensional 



5 

velocity dispersion ~ \/3 x 400 = 700 km s~^ between 
z = 2 and z = 4. 

In the PS formalism, the comoving number density of 
halos with velocity dispersion cr„ can be calculated by 

d\na{R,z) Sjz) / (5g(z) \ dR 

^ dliiR a{R,z)'''''^\ cj^{R,z)) R ^ ' 

As in Eq.(4), 6c{z) in Eq.(ll) is the critical overdensity 
for collapsed at redshift z. The meaning of R here is the 
same as in §2. The relationship between a^, and R is given 
by (Narayan k White 1988) 

o-„ = c„V3HoR{l + zfl'^ (12) 
for Einstein-de-Sitter universe, and 

CT„ = c^sflVL^'^H^Ril + z)i/2_ (13) 

for open universe or the LCDM universe. The coefficient 
Co- has been determined by comparing Eqs. (11) and (12) 
with N-body simulation (Jing & Fang 1994). It found c„ = 
1.1 — 1.3. For most calculations, Co- = 1.2 is a preferred 
value. Hence, one can safely use Eqs. (11)- (13) even when 
the 3-D velocity dispersion is as large as about 800 km 
s-\ 

The total number density of the collapsed halos with 
the velocity dispersion greater than a certain value, say 
cr''™, is 

N{>a''"',z)^ n{a^,z)day. (14) 

Using Eq.(14), we found that the SCDM and LCDM 
models can produce sufficient number of a'*™ = 700 km 
s~^ halos to fit with the number density of QSOs. But the 
number of such halos in the CHDM model is too small. 
To have enough number of halos, we should use cr'"" = 
300 km s"^ or less in the CHDM model. 

The mass of the (Tt,-sclcctcd halos can be determined 
from Eqs. (12) or (13). One can then calculate the two- 
point correlation functions of them by Eq.(lO). Here we 
get almost the same results as before because ay is one- 
to-one related to the mass or the radius of the top-hat 
windows. Similar to Fig. 2c, the cr'*™ = 300 km s~^ halos 
still lack of correlation power on all scales larger than 2 
h-i Mpc in the CHDM model. 

This result should be expected. An identification of 
objects from a density field is a sampling. As one knows, 
a sampled field will not be different from the original field 
on scales much larger than the characteristic scale, Rc, of 
the sampling (Vanmarkc 1983). This means that, in prin- 
cipally, no biasing on scales larger than the characteristic 
scale can be introduced by the identification. For our ques- 
tion here, it is not easy to given the characteristic scale Rc 
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of the identification. However, it is probably reasonable to 

choose Rc = [1 + ^{R,z)]'^/^R, where ^(i?, z) is the aver- 
age of ^(r, z) inside R (Hamilton et al. 1991). So the scale 
Rc is actually the coherence length of the density field 
smoothed by window R. The distribution of halos should 
not be biased from mass distribution on scales larger than 
Rf.. The correlation function can effectively be amplified 
by linear geometrical biasing only for large halos (Kaiser, 
1984). Figure 2 show a break in the correlation functions 
at rfe, which is roughly equal to Rc- When r > rb, the cor- 
relation functions of the collapsed halos should approach 
to the mass correlations, while when r < rb, the biasing 
leads to larger correlation functions. 

All identifications based on gravitational parameters, 
such as circular velocity, virial temperature etc. are essen- 
tially equal to mass identifications, because for a gravita- 
tional collapsed systems these parameters are one-to-one 
related. Therefore, they can be expressed as an equivalence 
to M. Biases in these parameters cannot give significant 
different results from those using cr^, or M. 

Obviously, any non-gravitational identifications are 
not constrained by the gravitational characteristic scale. 
For instance, if gas processes in the formation of QSOs 
play the role of biasing, the characteristic scale will not 
equal to Rc- However, the characteristic scale of gas pro- 
cesses should be much less than Rc, because the velocity 
of gaseous component is too small to segregate the QSO 
halos on scales of 5 h~^ Mpc. 

4. Conclusions 

We showed that the abundancc-plus-corrclation test can 
serve as a promising discrimination among models of cos- 
mic structure formation using the high redshift {z > 2) 
QSO number density and their clustering on large scales 
(10 h-i Mpc). We found that the SCDM and LCDM 
model are consistent with the observations of abundance 
and two-point correlation, while the CHDM models is dif- 
ficult to fit with the observed numbers. 

The contributions of variously biasing mechanisms 
have been investigated. A bias, in fact, is a model of the en- 
vironment suitable to form objects, or a phenomenological 
relationship between the cosmic density field and QSOs. 
We have studied possible biasing processes, including both 
gravitational and non-gravitational origins. None of them 
seems to give large characteristic scale and clustering am- 
plitude required to make the CHDM model success. 

One should also consider the possibility that each halo 
may host more than one QSOs. In this case, the correlation 
function can be stronger than that of one-QSO/one-halo 
model. In order to fit with observed abundance of QSOs 
in the CHDM model, one can assumed that each halo 
of M = IO^'^Mq host more than 10 QSOs. In this case, 
halos of M = 10"Mq in the CHDM can also fit with 
the observed two-point correlation function of QSOs (Fig. 
2c). However, this apparently leads to the paradox that 



we did not found so many bright QSOs concentrating in 

galaxy groups. A possible way to explain it is that the 
lifetime of QSOs is just longer than the cosmic time at 
z = 3 but shorter than the cosmic time at the present. 
(If QSOs have a life time shorter than the cosmic time at 
z = 3, we require even more low mass halos, this makes 
the situation even worse.) However, the difference between 
the cosmic time at z = 3 and z = is only a factor of 8 
{t (X {1 + 2)"^/^ in the = 1 universe). Therefore, we may 
predict that on average, each of 8 galaxy group will have a 
bright QSO. This is not true. Otherwise, we should require 
a very particular mechanism to make the QSO formation 
at = 3. 

It is theoretically possible to explain any correlations if 
we are allowed to introduce unknown inhomogeneity into 
the density field, and assume that the correlation of QSOs 
absorbers is given by these inhomogeneities. However, to 
plan these inhomogeneities is equal to put desired struc- 
tures in the initial perturbations. Such models will, how- 
ever, no longer be the SCDM, LCDM or CHDM models, 
which are based on initial fluctuations produced at the 
inflationary era. However, even in this case the CHDM 
will still be difficult to produce enough halos with reason- 
able velocity dispersions to fit with the abundance of high 
redshift QSOs. 
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